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Abstract 
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I. INTRODUCTION 



Although charge independence and charge symmetry are respected approximately by the 
strong interactions, these symmetries are both broken by the electromagnetic interaction 
and the mass difference of the up and down quarks. The symmetry-breaking remains of 
considerable interest both experimentally and theoretically because it can be studied by 
perturbation theory. At low energies, the large scattering length of two nucleons in the 
^5*0 state makes it particularly suitable for investigations of both charge independence and 
charge symmetry |]T]|. At medium energies, polarization experiments in n-p scattering have 
been of interest for studies of charge symmetry In nuclei, the mass difference between 
light-to-medium conjugate nuclear pairs (with a neutron replaced by a proton) also shows 
charge symmetry-breaking beyond that due to Coulomb forces; this is commonly referred 
to as the Nolen-Schiffer or Okamoto-Nolen-Schiffer anomaly @,|^. 

Theoretically, charge- dependent forces have been cast into 4 classes 0. It is the class 
IV forces, proportional to {r^^^ — T2^^) and (ri x T2)^^\ where r/ is the jth component of the 
Pauli isospin matrix for particle i, which produce charge-symmetry-breaking effects in the n- 
p system. In meson-theoretic investigations , these forces are due to (in addition to photon 
exchange) pion exchange, mixed p — uj and tt — t] — t]' exchanges, combined -k — ^ and other 
complex exchanges. The experimental asymetries measured at 477MeV and 183MeV 
are particularly sensitive to pion-exchange and p — uj mixed exchange, respectively. Indeed, 
both experiments agree with theoretical predictions based on meson theory 0. 

At 183MeV, the dominant contribution to the polarization asymmetry in the n-p elastic 
scattering cross section is found to arise from p — uj mixing in the exchanged mesons, as 
shown in fig.|l|. However, it has been pointed out by Goldman, Henderson and Thomas |§] 
that the fit makes use of the on-mass-shell mixing of the p and cu, as measured in e^e~ 
collisions, whereas it is space-like values of the 4-momentum transfer {(f < 0) that are 
required for the fit to the n-p scattering data. Furthermore, they point out that in a model 
which uses free quark-antiquark intermediate states, there is considerable variation of the 
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rho-omega mixing parameter with g^; indeed, it changes sign at ~ — (400MeV)^ and is 
small for negative (space- like) values of q^. This casts doubt on the theoretical calculations 
1^ . More recently, Piekarewicz and Williams have used an intermediate NN to calculate 
the mixing and obtain qualitatively similar conclusions as Goldman et al. 

The Nolen-Schiffer anomaly has also seen a revival of interest. In a meson-theoretic 
investigation, the type III forces proportional to {rf^ + Tg^^), in particular the p — u mixing 
contribution, is claimed to be an essential ingredient to explain the anomaly [l^. However, 



the on-mass-shell value of the p — uj mixing is used also in this analysis. The partial restora- 
tion of chiral symmetry in nuclei and the associated change of the n-p mass difference in 
medium is another explanation of the anomaly without recourse to the rho-omega mixing 
TT| . The effective theory of QCD |Tl| and QCD sum rules in a nuclear medium |T^ 



seem 



to support this explanation. Further extension of the approach in ref. [T^ is also made by 
taking into account the mixing in the vector channel [1^. Thus it is important to find out 
the relative importance of the p — uj mixing among other explanations by analyzing the 
variation of the mixing. 

Of the various non-perturbative approaches in QCD, QCD sum rules have been found to 
be a powerful tool for analyzing the properties and decays of hadrons 0,|TB[ . They can give 



a plausible explanation of the Nolen-Schiffer anomaly |]T2| as we mentioned above. One of 
the earliest applications of the sum rules was actually to a calculation of the on-mass shell 
p — ijj mixing . 



In this paper, we follow the path forged by Shifman, Vainshtein and Zakharov (SVZ) 
||16|| , but we extend it to space-like values of which is the relevant region of the charge- 
dependent nuclear forces. Furthermore, we re-evaluate the various semi-phenomenological 
constants required by the sum rules. Whereas SVZ used a value of (5 (see below) obtained 
from a "speculation" based on "the experience with sum rules" , or as a free parameter, we 
obtain (5 from the sum rules themselves. In section 2, we will summarize the effect of the 
p — UJ mixing on the nuclear force and give an essential idea about the importance of the 
g^-dependence of the mixing. In section 3, we will determine the parameters which control 



the variation from QCD sum rules. The combined use of the Borel sum rules and the 
finite energy sum rules together with the current knowledge of the quark mass difference and 
the quark condensates give us an unambiguous determination of these parameters without 
further assumptions. We find that some of the outputs are different in magnitude and sign 
from the SVZ results. In section 4, the dependence of the p-uo mixing and also its effect 
on the nuclear force are discussed by using the parameters determined in section 3. The 
effective mixing angle 9{q^) changes sign at q^ ~ OAm^p and gets negative in the space like 
region, which is qualitatively consistent with the previous analysis ; however the effect is 
stronger in our case. In coordinate space, this corresponds to a node of the charge-dependent 
potential around r ~ 0.9fm which suppresses the effect of the p — uo mixing in the nuclear 
force. Section 5 is devoted to a summary and concluding remarks. 



11. p-u MIXING IN THE SPACE LIKE REGION 

The Feynman graph for the class IV forces originating from p—uj mixing is shown in fig.|l|. 
The cross in fig.|l| denotes the p — uj mixing 9{q^) defined in terms of the mixed propagator 

^7M^) = ^ I c?W''-(Tp^(x)c.,(0))o, (2.1) 



\9^lu r, j - 



where we have neglected the width of the resonances for simplicity. 

At q^ = m^, 6{q^) reduces to the mixing matrix (p | Hsb \ 'jj) measured in the e^e~ — *• 
TT+TT" data 

0{q' = ml) = {p I Hsb I oo) 

= (-4520 ± 600) MeV^. (2.2) 

Here Hsb is a part of the Hamiltonian which breaks isospin symmetry and we have used the 
covariant definition of the mixing matrix (p^ | Hsb \ ^u) = —{Ofiu — (liilu/ q^){p \ Hsb \ ^) 



following ref. |T8 



HsB has both QED and QCD origins. One of the QED effects to eq.(2.1); from p — > 
7 — tu; can be exphcitly calculated and gives a small and positive contribution p!9[ 

= - 610 MeV^, (2.3) 

where e/gp and e/g^ denote,respectively, the electromagnetic coupling of p and u with 7. 
We have assumed (7^ ^ 3(7p,as obtained from the flavor SU(3) symmetry. Thus the main 
part of the observed mixing comes from the quark mass difference 

"^d - - 0.28(md + m„), (2.4) 

where (m^ + m„)/2 = (7 ± 2)MeV at a IGeV QCD scale ( [|0[ and ref.therein). Note also 
that, once one takes into account the 7 exchange between nucleons with proper nucleon 
electro-magnetic form factors (fig.0) , one has to subtract the p — 7 ^ effect from fig.[l| to 
avoid double counting. (The QED mixing is hidden in the form factors in fig.0.) 

For isospin mixing in the nuclear force, we are concerned with space-like momenta < 0, 
which is far from the on shell point = m'^p^. Thus it is crucial to study the variation of 
6{q'^) to determine whether one can naively utilize the on-shell mixing value, eq.(2.2), in the 
space like region |Q] . In the next section, we will analyze the variation of p — u mixing on 
the basis of the QCD sum rules in which the physical quantities are related to the condensates 
in the QCD vacuum. Our approach has less ambiguities compared to other approaches based 
on phenomenological models Before entering into details, let us discuss the primary 

origin of the variation to demonstrate the essential idea in a model-independent way. 

We start with the unsubtracted dispersion relation for 11^'^ (g^): 

^^K%^') = - '-^ds. 2.5 

Ti Jo s — 

If g^ is not too large and in the space-like region, the main contribution to Rell^'^ comes 
from the lowest poles in the imaginary part of f^^, i.e. the p and u resonances. It is 
important to note that for the mixing, H^"^, unlike the case for HJ^^, the integral of the 
continuum contribution such as that from an NN loop is finite M. If photon exchange is 



included in the loop, this is no longer the case; however, we find this contribution is to be 
small numerically (~ 2%), so that it can be neglected. Thus we do not need any subtraction 
in the right hand side of eq.(2.4) as far as the strong interaction is concerned. In the pole 
approximation, lmll'^{s) is easily evaluated by saturating the intermediate states by p- and 
a;-mesons: 

^lmn^^(s) = n{0\p^\p){p\u,\0)6{s - mj) + 7r{0\p^\u) {u;\u,\0)6{s - ml) 

= 7TFp5is - ml) - TcFJis - ml). (2.6) 

Since mp ^ m^, the residues at the p and cj-poles do not have to be equal 

Fp-F^ = 0{md-m^). (2.7) 

Substituting (2.5) into (2.4), one obtains 

'Ren-(g^) = ^1 

mp 



3 — g2 — ml ' 



5m^{Fp + F^)/2 - (g2 _ m^)(Fp - F^) 



(2.8) 



[q^ — m^p){q^ — ml) ' 

where = (m^ + ml)/ 2 and 5m^ = ml — = 0{mfi — m„). 
Thus we finally get 

%') = ^K)[1 + A(4-1)]' (2.9) 

with d{m^) = 5m^{Fp + F^)/2 and A = -{Fp - F^)m^/e{m^). The first term of eq.(2.8) is 
nothing but the usual on-shell mixing matrix {p\Hsb\^^) , while the second term gives a 
variation due to Fp^ F^^. Since both the first and second terms are proportional to m^ — mu 
and are of the same order, there is no a priori reason to neglect the second term. However, 
in most phenomenological applications [|l],^|lO| of p — uj mixing to the charge symmetry 
breaking nuclear force, Fp = F^^ has been assumed without any justification. 

The variation of 6{q'^) is an inevitable consequence of the isospin symmetry breaking 
and we have not used any specific models up to this stage. However, the sign and the 



6 



magnitude of Fp — F^^ (or A ) in eq.(2.8) depends on the QCD dynamics. In QCD, we have 
to start with the quark correlator instead of the hadronic correlator 11^'^: 

n^.(g') = ^ J d'xe''^^{TJP{x)J::mo, (2.10) 

where 

= (m7^m - d-f^d)/2 j; = {u-fpu + d-f^d)/6. (2.11) 

These vector currents couple to p and u as well as to the higher resonances {p',u!',- ■ ■) 
and the continuum. Thus the imaginary part of U^i, is written as the sum of all the mixed 
hadron correlators 

lmU,,is) = AoImn^^(.) + Ailmn^';r'(.) + ■ • ■ , (2.12) 

where An{n = 0,1,2, ■■ ■) denotes the overlap of the quark current with the physical hadrons. 
For example 

2 2 

mm, 

Ao = (2.13) 

9p9u; 

where gp^^ is defined by (0| J^''^|p(u;)) = {m^p^^/ gp^^)tp with g^ ~ S^fp as given by flavor SU(3) 
and gl/A-n ~ 2.4. 

By using the dispersion relation for 11^^ and the OPE for Ren^^(g^) at — ^ — oo, one 
can extract the resonance parameters such as the mixing matrices for p — uj and for higher 
resonances in the r.h.s. of eq.(2.10). The first attempt at such an analysis was made by 
Shifman, Vainshtein and Zakharov |16|. However, one of the mixing parameters (3, which 
is related to A as A = [3 + 1, was not determined well because of the poor e"'"e~ data 
available at that time and of the large uncertainty of the current quark masses. In the next 
section, armed with the updated information of these parameters (eqs. (2.2) and (2.3)), we 
will determine the crucial quantity (3 within the framework of the QCD sum rules. The 
combined use of the finite energy sum rule and the Borel sum rule is a key ingredient there. 
From the analysis, we find that A is indeed positive and, including all uncertainties in the 
parameters, is found to lie in an interval: 
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A = /5 + 1 e [1.43, 1.85] (2.14) 



(5 is close to, but generally larger than the value assumed in ref. 16 , 0.5. It follows that 



^(g^) has significant variation in the space like region; indeed, it vanishes at 

q^ = m^^^_ (2.15) 
A 

Thus the nuclear force is significantly affected by the variation of 0{q^). This can be 
illustrated by examining the coordinate space charge-symmetry-breaking nuclear force in 
the static limit which is related to the Fourier transform of eq. (2.1): 

= -M^^fi _ ^]e-^ (2.16) 
47r 2m mr 

where gp{uj)N denotes the vector coupling constant of the nucleon and vector mesons. (We 
did not consider the form factor for the coupling constants, because we are interested in the 
longer range part of V^.) When (3 = 0, eq. (2.11) reduces to the potential given by Coon 



and Barret |]T^. The potential is strongly attenuated by the presence of (3 and becomes zero 
at 

2A 

r = 0.9 fm , (2.17) 

m 

which is the region of interest for symmetry-breaking effect. 

III. p-uj MIXING IN THE QCD SUM RULES 

A. Operator Product Expansion (OPE) 

For the QCD sum rule determination of the p—uj mixing, one should start with the current 
correlator |[T^ given in eqs. ( p.lO| , p.llD . Because the currents and are conserved, II^i, 



has to be of transversal structure: 

n^.(g') = (»-^?M-g')n(g') (3.1) 
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Following ref. [|T6|, we will analyze sum rules obtained by the OPE for n(g^). In the OPE 



series, we keep power corrections up to order 6. Furthermore, in each term we take con- 
tributions which are of 0{mg) {q = u or d) or of 0{a), and any terms 0{m?^ , 0{a'^) and 
0{mqa) are neglected. We will check the validity of these assumptions later on. Therefore 
we are faced with the diagrams in fig.^ which contribute to Il{Q^) (Q^ = —q^) as follows 

(1) Free Part (fig. 3a): The leading order term proportional to InQ^ vanishes, while the next 
term 

_ 3 mj-ml 

- 2^^2Q^' ^^-^^ 
is of 0{m'^) and will be neglected. For the same reason, the one-gluon exchange diagram 
(fig. 3b) is of higher order in rriq. 

(2) Irreducible 1 — 7 exchange (fig.3c): Because of the different coupling of the 7 to the u 
and d quarks, the one-7 exchange in fig. 3c contributes as 

U(2) = lnQ2. (3.3) 

167r3l2 ^ ^ ^ 

(3) Reducible 1 — 7 exchange (fig. 3d): As has been proposed in ref. it is convenient to 



separate the reducible 7 exchange diagram from the whole analysis from the very begin- 
ning. This can be done by neglecting fig. 3d in the OPE side and subtracting p ^ 7 — > u; 
contribution (|2.3| ) from the phenomenological side. 

{4)Gluon Condensate as{G'^)o (fig.3e) : It vanishes for the same reason as (1). 
(5) Quark condensate mq{qq)o (fig.3f): 

2 _ 

n(5) = Y2Q4f"^''^""'^° ~ md{dd)o] . (3.4) 



(6) Mixed Condensate (fig.3g): As has been shown in refs. p!6|j21[] , even in case of only one 
quark flavor, these condensates contribute at least in (9(mg) to the vector current and are 
therefore not present in our analysis. 

(7) 4-Quark condensate (fig.3h-k): There are contributions from gluon as well as 7 exchanges. 
Typical diagrams are shown in figs.3h-k. 
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(a) (fig.Sh) 

n(7a)(Q') = 

(b) (fig.3i) 

n(7fc)(Q') 

(c) (fig.sj) 

n(7c)(Q') 

(d) (fig.3k) 

n(7d)(Q') 

B. Symmetry-breaking parameters and vacuum condensates 

There are essentially two quantities which control the isospin-symmetry breaking in OPE: 
fndlmu and {dd)f)/ {uu)q. Although they are related in principle in QCD, they appear as 
independent parameters in the QCD sum rules. As we will see below, the dimension 4 
operator in the OPE is controlled by the quark mass ratio while the dimension 6 operator 
depends on the ratio of the condensates. Because the dimension 4 operator is the dominant 
term in the OPE, the p — uo mixing is essentially determined by the quark mass ratio and 
depends weakly on the ratio of the condensates. 

The quark mass ratio is determined by the analysis of the mass splitting of mesons and 
barons using the chiral perturbation theory ( and ref .therein) : 

— ^=0.28 ±0.03 . 3.9) 

nid + ruu 

The isospin breaking of the quark condensate has been analyzed by several methods; the 
chiral perturbation theory the QCD sum rules for scalar and pseudoscalar mesons 



(3.5) 



Ana, 



9 ■ 12Q' 



(3.6) 



Sna 4 2 1 ^ y ,n2\ 

Y^(g(w7a75M) --^{dlal^d) )o 



(3.7) 



Wna ,4 2 ,\2\ 



(3.8) 
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5[, effective models of QCD incorporating the dynamical breaking of chiral symmetry 
26| , |27[] and the QCD sum rules for baryons p8 |. 



7 



{dd)o 
{uu)o 



-(6 ~ 10)10-3 [22] 
(10 ±3)10-3 [gl 



(3.10) 



MM 



-(7-9)10-3 
-(2 ± 1)10-3 

Here we have assumed {ss)q/ {uu)q — 1 = —(15 ~ 30)% to get the first number in the above. 
The last number is obtained solely by using the baryon mass splittings in the QCD sum rule 
and is quite different from the other determinations. 

Besides the symmetry-breaking parameters, we need to fix the value of the four-quark 
condensate. The vacuum saturation hypothesis gives a rough estimate of its magnitude; it 



gives a satisfactory description of the p-meson properties |M. However, there are recent 
arguments which favor a value larger by a factor 2 or more than that adopted in (see 



e.g. |]T5|,^ and the references therein). Since this condensate is not determined well in the 
sum rule, we will take the following two typical numbers and carry out the Borel analysis 
using them: 



2 1.81-10 '^GeV^ vacuum saturation [14] 
3.81 • 10-^GeV^ [H 



(3.11) 



The definition of 7 in eq.( |3.10| ) and the vacuum saturation hypothesis allow us to rewrite 
the isospin breaking in the four-quark condensate as 



{uuy, - {dd)i = -2^{qq)i 



(3.12) 



a{uu)1 = a{dd)l = a{qq)l 



(3.13) 



Using the above symmetry-breaking parameters and the vacuum condensates together 



with the Gell-Mann-Oakes-Renner relation ||30| 



{rriu + md){uu + dd)o = 



(3.14) 
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eq. (|3.4] ) and the sum of ( |3.5| )-(p78D are written as 



n(5)(g' 



1 md-mu^,2 



ijid + rriu 



1 + 



md + rriu 7 
mrf - m„ 2 + 7 



(3.15) 



and 



n(7)(Q' 



224 2 



81 i2g 

Summing up all contributions, we obtain 



-7 + 



a 



(3.16) 



12n(Q^) = -colng^ + |L + ^ + 2^ 



(3.17) 



with 



a 



Co 



Cl 



IGvr^ 



md-mu 2 2 
C2 = : 



md + rriu 
224 



1 + 



C3 



^1 



-vr 



as(gg)o 



md + rriu 7 
TTid - m„ 2 + 7 
a 

-7 + 



8Q;s(Ai^) 



(3.18) 
(3.19) 
(3.20) 

(3.21) 



For the renormalization point /i^, we take a typical scale of the Borel mass IGeV^ so that 



as = a,(lGeV^) ~ 0.5 . 



(3.22) 



It is obvious from eqs. ( p.20| , p.21| ) that the md/rriu controls the magnitude of the dim. 4 
matrix element and {dd)^/ {'Uu)o controls the dim. 6 matrix element. 

The aim of our sum rule analysis is to determine the parameter (3 (or A); we will use 
the experimental p — 00 mixing at = as an input. For later use, let's define here the 
dimensionless mixing of hadronic origin as 

-12 ^fm^) ^ 



^11.3 ±1.3) -10" 



(3.23) 



To take into account the large uncertainties of 7 and as(^5')o, perform our analysis 
for the 4 different sets of input parameters specified in table| with the corresponding OPE 
coefficients q (i = . . . 3) shown in table|I|. Borel stability analyses will be made for these 



four different sets in section III E 
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C. Spectral Function and Sum Rule 



The general expression for the spectral function, i.e. the imaginary part of n^,y(g^) is 
given in eq.( 2.12 ). As has been pointed out in refs. [rB,r3|, it is absolutely necessary to keep 



at least one higher resonance (p', uj') in ( p^.l2[ ) in order to obtain stability in the QCD sum 
rule. We will treat the related mass m'^ = |(m^, +m^,) as well as the corresponding coupling 
constants Qpi and g^^l as effective parameters to be determined from the Borel stability 
analysis. Furthermore, in order to account for the electromagnetic logarithm in eq.( p.3|) in 
the space-like region, we have to introduce a continuum threshold sq. The corresponding 
part in the spectral function is 'suppressed' by the electromagnetic coupling a and indeed it 
will turn out that the influence of this term on the sum rule analysis is very small. Assuming 
sharp resonances for p^uj^p' and uj' as well as a step function for the continuum one can write. 



analogous to eq.( |2.6| ): 



12Imn(s) = 7r/,5(s - ml) - nfj{s - ml) 



+ 7r/p'5(s - mj) - nf^>6{s - m^,) 



The general form of the Borel sum rule reads 



/ dse-^/*' Imn(s) = nB(M2), (3.25) 
Jo 



7rM2 Jo 

where the Borel transform 11^ (M^) is defined by: 

n.(M^)^ lim ^Q-(|,)"nW^). (3.26) 

Q^, n — ^ oo 

M2 = gVn fixed 
Let us define the following parameters 

5m'^ fp + fu^ 



2 



_ Sm'^ fp' + L' 
^ ~ m'' 2 ' ^^-^^^ 
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5m 



2 /p "I" A; 



5m 



,2 fp' + fuj 



(3.28) 



^ and ^' are directly related to the on-mass-sliell mixing in the p — u and p' — uo' channels, 
respectively. (3 and (3' are related to the isospin-symmetry breaking at the resonance poles 
and control the dependence of the p — uj and p' — uj' mixing. The relation of (3 and 
A defined in section 2 is easily obtained by comparing ( |2.12|) , (|2.13|) , (|2.6|) with and 
(H): 

A = /3 + l . (3.29) 

By using the above definitions, one obtains from eq. (|3.17 ) up to 0{5m'^) and 0{5m''^) 
the Borel sum rule (BSR): 



(3.30) 



where C0-C3 are in units of 10^"^. 

Starting from this expression, one can carry out the Borel analysis (stability analysis 
with respect to M^) or analysis of the finite energy sum rule (expansion in l/M^ on both 
sides). Let us summarize our strategy before entering into the details. On the right hand 
side of eq.( p.30|) , we have 7 resonance parameters: /3, /?', ^, m, m' and sq. Among them, 
m is the average mass of p and uj, and ^ is directly related to the measured p — uj mixing, 
so they are known inputs, m' should be chosen to be a number close to mpi = 1390MeV or 
m^/ = 1450MeV; thus it is not arbitrary. The sum rule depends on sq very weakly because 
of the suppression by a; thus the choice of Sq does not materially affect the result. Hence 
we are left with essentially three unknown parameters (3, (3' and ^' to be determined by the 
sum rule. 

In the finite energy sum rule, by expanding eq.( |3.30| ) in powers of l/M^, we get three 
independent equations. They are enough to solve for the three unknowns. One can also 
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change m' and sq slightly to check the sensitivity of the result to that change. One of the 
problems of the FESR when compared to the Borel sum rule is that it is generally sensitive 
to the parameters of the higher resonances. Therefore, after getting a rough idea for the 
magnitude and sign of the three unknowns, we will move on to the Borel sum rule. 

By taking a derivative of eq.( |3.3(JD with respect to M^, one generates another independent 
sum rule. To carry out the stability analysis for each unknown parameter (/?, (3' and ^'), we 
need one more constraint, and we take a first sum rule of the FESR to supply this constraint. 
The first sum rule is known to be a local duality relation and is least sensitive to the higher 
resonances. By this procedure, we can draw independent Borel curves for each unknown 
and carry out the stability analysis by varying m' and Sq, i.e. ml and Sq are chosen such that 
/?, j3' and ^' are least sensitive to the variation of the Borel mass M^. In this way, we can 
completely determine the resonance parameters within the QCD sum rules. Our procedure 
is different from the original SVZ analysis where (3 was assumed to be 1/2 as an input 
and also poor experimental data of was used. 



D. Finite Energy Sum Rule (FESR) 

The FESR can be formally obtained by expanding the BSR eq. (|3.30|) in powers of 
and comparing the coefficients of each power on both sides. In doing so one obtains three 
FESR: 

- /5e - ^ - co4 = (3.31) 
« + ^^ + ^'^+?'/3'^ + lc„(il)=^c. (3.32) 

2 mP 2 rrr o \m 



< - -.iP - e— - -J'P'-^ - £^ ) = C3 (3.33) 



Using various values for m'^ and sq as input parameters we can therefore solve the equations 



( p.31| )-( p3^ ) which gives us the 3 unknown parameters (3, ^' and jS'. TablePTT] shows their 



values for the OPE-set I. As we can see, the result is sensitive to the choice of the resonance 
mass m'^, whereas the influence of the electro magnetically suppressed threshold Sq is small. 
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It is important to notice that [3^ and j3'^' have to be opposite in sign due to the first sum 
rule of FESR ( |3.31|) or equivalently the local duality relation, 

roo 

/ dslmliis) = 0(a) ^ , (3.34) 
Jo 

This is the reason why one should keep the resonance explicitly |T^, otherwise the term (3^ 
cannot be canceled. 



E. Borel Analysis 



Finally we perform the Borel stability analysis of ( |3.30|) . Again we take for ^ the exper- 
imental input ( p.23|) . m'^ and sq are used as variable parameters which are fixed in order 
to obtain maximal stability of the Borel curves. In order to choose a stability window for 
M^, we demand that, for each parameter set (I-IV), the contribution of the 6th order power 
correction in the OPE is less than 25% of the 4th order one. This gives a lower limit M^j„ 
for the Borel mass M^. 

As we mentioned before, we adopt the first sum rule of the FESR ( p.31| ) or equivalently 
the local duality relation ( p.34|) as an additional constraint. Also, we generate another sum 
rule from (|3.30|) by operating with jji both sides. In doing so one can 

uniquely calculate the 3 unknown parameters (3, ^' and (3' as function of M^. 

1 Z2(M2) + Zi(M2) (l-^) 



(3{M^) 
(3\M^) 



eiyi(M2)(^-i)-i^2(M2 



'M 



2\ 2 



m 



l2 



So m 
m m' 



(3.35) 

(3.36) 

(3.37) 
(3.38) 



with: 



Zi(M2 



Co 



+ Ci 



m 



(1- 
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m 



2\ 2 



2 



+ C2 



m 



2 \ 2 



+ C3 



m 
M2 



2\ 3 



m 
JP 

1 
M4 

1 

M4 
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9(1/M2) V'"' '"^i^'"' 



a(l/M2 
5 



(3.39) 



The general procedure for determining m'^ and Sq is now the following: For each of the 
parameter sets I-IV from table|l| we first take the corresponding central value (i.e. without 
error bars) of the OPE coefficients q (i = 1 . . . 3) from table|l| and calculate the Borel curves 
/3(M2), ^'(M2) and (3'{M^) due to eqs. (|06|) , (|07D and (|08|) , respectively, for various 
values of m'^ and sq. The final result is the one with the least dependence on within 
the Borel window > M^i-r^ with Mmin^ determined above. The results for /5 and ^' 
using parameter set I are shown in figs. ^ and ^, respectively. One recognizes a wide stable 
region, which even can be extended up to ^ lOGeV^. The maximal stability occurs 
for m'^ = 1.6GeV^, where Sq can be increased up to about 2.4GeV^ without any significant 
change. This insensitivity of the result to the variation of the continuum threshold are 
found for all other parameter sets as well and refiects the fact, that the dependence on sq 
is suppressed by the electromagnetic coupling a. With m'^ and Sq obtained in this way, we 
now consider the change in the Borel curves for (3, ^' and (3' for the upper and lower limits of 
the Ci from table|| and extract from this the upper and lower limits for (3, ^' and For f3 
(set I) the limits are shown in fig.^ The final values are given in table^. It turns out that 
the variation of /3, ^' and (3' is completely determined by the variation of the coefficient C2, 
which itself is based on the uncertainty of the quark mass difference = 0.28 ± 0.03. 

The uncertainty of 7, which creeps into the coefficient C3 of the dimension 6 condensate 
can be neglected. It should be noted that a full readjustment of m'^ and sq for the upper 
and lower limits of the q's from table|I| which are obtained by performing a new stability 
analysis for each case would result in much bigger error bars for (3, ^' and (3'. 

Similar considerations hold for the parameter sets III and IV. In case of set II we find 
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that the Borel curves are rather unstable in the window > Mmin^ = l.SGeV^ (figs.0 
and H), especially the one for and therefore this uncertainty, which is bigger than the 
variation due to the Cj's, has to be included in the corresponding error bars of table|V[ 

Finally we have convinced ourselves that in none of these cases does the result change 
noticeably if one either switches off the electromagnetic interaction completely, i.e. puts a 
equal to zero, or includes the 0{mq^) diagram (fig. 3a) with the corresponding coefficient ci 
in the Borel analysis; this analysis justifies the neglect of these terms in higher orders. 

We also want to stress that the value of (3 obtained by this method differs significantly 



from that used in ref. (3 = 0.5, which was based on a SU{3) symmetry argument and 
served as input parameter. Also the individual signs of ^' and /?' come out different than 



m 



16 1 , which arises from the fact, that in addition to /3, a different experimental value of ^ 



was used there. 



IV. 6i(q2) and the n -n potential 

In the previous section, we have determined the parameter (3 (=A — 1) which controls 
the g^-dependence of the p — uj mixing parameter d{q^). The result summarized in table(rV| 
together with the relation between A and 9{q^) (eq.( p.9| )) immediately leads to fig.^. The 
region inside the shaded band of fig.^ is favored by our QCD sum rule analyses. 

An important observation is that A is always greater than 1 as long as (3 is positive. This 
feature leads to the conclusion that 9{q^) already changes sign for positive and is always 
negative for < 0, as can be seen in fig.|^. Although the considerable variation of 9{q'^) 
obtained here is qualitatively similar to previous results based on different models, our d{q^) 
has a stronger dependence than that of others. 

In the quark-loop model {p ^ qq ^ uj) with momentum cutoff, Goldman et al. [§, 
^(g^) changes sign at a space like momentum (g^ ~ — m^/2), therefore the variation of ^(g^) 
is more moderate than ours. The nucleon-loop model (p — > NN —>■ uj) with dimensional 
regularization predicts a sign change of ^(g^) at g^ = 0, which still corresponds to a more 
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moderate variation than our result. 

One should notice here that there is a crucial assumption with no theoretical justification 
in both quark-loop and nucleon-loop models: the effect of the isospin breaking other than the 
QED effect is solely attributed to the mass difference between u and d constituent quarks (in 
the quark-loop model) or to the mass difference between the proton and the neutron (in the 
nucleon-loop model). There is no a priori reason, however, to neglect the isospin breaking 
in the coupling constants of the vector mesons with the constituent quarks or the nucleons, 
which generates an extra effect to the p-u mixing of 0{md — friu); nor is there justification 
for neglecting nuclear resonances, e.g. A, N*. Thus it is fair to say that refs. take 
into account only a part of the total contribution due to the current-quark mass difference. 
We do not suffer from this deficiency, since all the isospin-breaking effects are automatically 
taken into account in the OPE of the current correlator. 

Let us also comment on the unsubtracted dispersion relation for 11^^ (g^) which we have 
used to extract 9{q'^) from Imn^^j(g^). As can be seen from the OPE in section 3, n^^(g^) 
approaches — as — > — oo if we switch off the small QED effect. Therefore, the 
dispersion integral is convergent and there is no need of a subtraction to extract Ren^^(g^) 
for the entire region from the spectral function obtained from the QCD sum rule. This is 
the reason why we could predict the g^ dependence of the mixing. If a subtraction were nec- 
essary, then extra experimental information would have been required to fix the subtraction 
constant. [] 

The g^ variation of ^(g^) is of primary interest for the isospin-breaking N — N interaction 
shown in fig.||. The N — N potential, including the effect of 6{q'^), is given by eq.(2.16). 
To compare the potential V^^{r) with and without the g^-dependence, we have shown it 
in fig.|l3 for X = P + 1 (the dotted line) and A = the (the full line), where we have used 



^When we determine the phenomenological parameters in the spectral function Imn^^(s) of the 
QCD sum rule, the subtraction is irrelevant since the Borel sum rule is used. 
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A = 0.78 from set I (c.f.table|rv|) as a typical example. The long range Yukawa potential due 
to the p — uj mixing (the first term in eq.(2.16)) is strongly suppressed by the dependence 
of the mixing (second term in eq.(2.16)). As a result, in the coordinate space the potential 
changes sign at r = 2A/m = 0.9fm. The large difference between the solid and the dashed 
curves of fig. at short distances is not really relevant, since the short range part of the 
potential V^j^{r) is screened by the strong central repulsion of the nuclear force and also is 
reduced by the the p — N — N and u — N — N form factors. 

Several remarks are in order here on the N — N potential due to the isospin breaking: 

1. In the conventional application of the p — u mixing to the N — N potential, only the 
mixing of the p — u propagators is taken into account, as is shown in fig.|l]. However, 
there is another source of isospin-breaking due to the isospin-breaking p{uj) — N — N 



coupling (c.f.figjTl|). The magnitude of this symmetry breaking is unknown. Further- 
more this effect does not contribute to the polarization difference in n — p scattering, 
but does contribute to the Nolen-Schiffer anomaly. 

2. The QCD sum rule can give us the mixing in the p — u channel and that in the p' — u' 
channel separately. Therefore we can study the effect of p' — u' mixing to the N — N 
force. From the QCD sum rule we have extracted ~ —2.5 ... — 3.5, so that the 
corresponding A' = + 1 is clearly negative and therefore the effect goes into the 
other direction, i.e. the p' — uj' mixing angle 0'{q^) increases in the space-like region. 
The effect of 9'{q'^) {p' — uj' mixing) to the nuclear force will appear only at short 
distances compared to 9{q'^) [p — uo mixing). Thus it is not so relevant to the low 
energy N — N scattering. Since we do not have any information of the p'{uj') — N — N 
coupling constants, we cannot make any quantitative estimate to the effect of 9'{q^) at 
the present stage. In the one boson exchange approach to the nuclear force the higher 
resonance contributions are not taken into account explicitly. 

3. In our approach we have taken a sharp resonance for the p meson spectral function 
and therefore completely neglected the width of the p meson. However we believe that 
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including the p meson width does not change our result significantly, because only 
the integral over the spectral function is important in the sum rule and furthermore 
the crucial quantity \ = (3 + 1 > 1 does not depend so much on the details of the 
calculation of (3. 



V. SUMMARY AND CONCLUSIONS 

QCD sum rules have been shown to be a powerful technique for studying meson and 
baryon properties. In this paper we have used these sum rules together with the usual 
unsubtracted dispersion relation to study p-uj mixing both on- and off- mass shell. We find 
a rapid variation of the mixing parameter d{q^) with . Thus, although we fit the value of 
6 a.t = nip^ ^ fn^'^, its sign has already changed at q^ = 0. Since it is the space-like value 
of that play a role in the contribution of p-u mixing to the nuclear force, our results have 
important implications for fitting experimental results. For example, the on-shell value of 
the p-u mixing parameter, ^(m^), has been used to fit the asymmetry observed in polarized 
n-p scattering [0-0,0] • Indeed, at ^ 200MeV, this (assumed on- mass-shell) mixing gives 
the dominant contribution to the observed asymmetry. Our result, like those obtained 
previously by other methods HJ^, thus place a large question mark on our understanding 
of the observed charge asymmetry. In addition p-uj mixing contributes to the difference of 
the p-p and n-n scattering lengths and to the energy differences of mirror nuclei. Indeed, 
there has been a revival of interest in the contribution of p-u mixing to the Nolen-Schiffer 
anomaly 0,0] and our result also puts some of these explanations [|l0|,p!3| into questions. 
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FIGURES 



FIG. 1. Class IV nuclear force from the p-uj mixing 
FIG. 2. 7-exchange between two nucleons with electromagnetic form factors 

FIG. 3. Diagrams in the OPE 
FIG. 4. Borel curve for /? with 3 different values of (m'^ , sq) (Parameter-Set I) 

FIG. 5. Borel curve for ^' with 3 different values of (m'^ , sq) (Parameter-Set I) 

FIG. 6. Variation of /3{M'^) due to the errorbars in table|| and tableHI where m'^ = 1.6GeV^ 
and So = l.SGeV^ are fixed and taken from fig.^ as the values with maximal Borel stability in the 
window M2 > M^m^ = l.OGeV^. 

FIG. 7. Borel curve for f3 with 3 different values of (m'^ , sq) (Parameter-Set II) 

FIG. 8. Borel curve for ^' with 3 different values of (m'^ , sq) (Parameter-Set II) 

FIG. 9. The off-sheh p-io mixing angle 6i(g2) with /3 = 0.78 

FIG. 10. The N — N Potential resulting from the p-uj mixing V^^{r) (in arbitrary units) 
calculated with A = 1.78 (dotted line) (which is the result for Parameter-Set I (c.f.tabl^)) in 
comparison with using the on shell value for (corresponding to A = 0) [|l^] (full line). The second 
figure magnifies the relevant region of r, where the sign change of the off shell curve takes place. 

FIG. 11. Isospin breaking nuclear force due to p or a; exchange with isospin breaking 
p{uj) — N — N vertex 
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TABLES 



TABLE I. The sets I-IV of input parameters for the OPE 




^•10^ 


mc^+mu 


7-10^ as( 


qq)l ■ lO^GeV^ 


Set I 
Set II 
Set III 
Set IV 


11.3 ± 1.3 
11.3 ± 1.3 
11.3 ± 1.3 
11.3 ± 1.3 


0.28 ± 0.03 
.28 ± 0.03 
0.28 ± 0.03 
0.28 ± 0.03 


-9.0 ± 3.0 
-9.0 ± 3.0 
-2.0 ± 1.0 
-2.0 ±1.0 


1.81 
3.81 
1.81 
3.81 


TABLE 11. 
tablel 


The OPE coefficients 


Cj (i = . . . 


3) calculated with the parameter 


sets I-IV from 




CO 


ci 






Set I 
Set II 
Set III 
Set IV 


0.147 
0.147 
0.147 
0.147 








2.66 ± 0.30 
2.66 ± 0.30 
2.70 ± 0.30 
2.70 ± 0.30 


- 0.83 ± 0.23 

- 1.75 ± 0.48 

- 0.29 ± 0.08 

- 0.61 ± 0.17 


TABLE III. j3,^' and j3' from the FESR for different values of m'^ and sq (Parameter-Setl) 


(m'^5o) [GeV2 


] 




i' 


P' 


(1.4 , 2.0) 
(1.6 , 2.0) 
(1.8 , 2.0) 
(2.0 , 2.0) 




0.91 
0.75 
0.64 
0.55 


1.03 
0.87 
0.74 
0.62 


-4.41 
-3.79 
-3.44 
-3.20 


(1.8 , 1.4) 
(1.8 , 2.0) 
(1.8 , 2.4) 




0.64 
0.64 
0.64 


0.73 
0.74 
0.73 


-3.38 
-3.44 
-3.52 
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TABLE IV. The final values of (3, ^' and P' as well as the resonance mass m' for the four 
different parameter sets of table| extracted from the Borel stability analysis 

/3 ^' /?' m'^ 

Set I 0.78 ± 0.07 0.98 ± 0.17 -3.52 ± 0.30 1.6 

Set II 0.65 ± 0.13 0.85 ± 0.46 -2.71 ± 0.60 2.0 

Set III 0.47 ± 0.04 0.56 ± 0.07 -2.61 ± 0.15 2.4 

Set IV 0.62 ± 0.06 0.76 ± 0.11 -3.06 ± 0.23 1.9 
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